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ABSTRACT 

It has been proposed that RX J1914. 4+2456 is a stellar binary system with an orbital 
period of 9.5 mins. As such it shares many similar properties with RX J0806. 3+1527 
(5.4 mins). However, while the X-ray spectrum of RX J0806. 3+1527 can be modelled 
using a simple absorbed blackbody, the X-ray spectrum of RX J1914. 4+2456 has 
proved difficult to fit using a physically plausible model. In this paper we re-examine 
the available X-ray spectra of RX J1914. 4+2456 taken using XMM-Newton. We find 
that the X-ray spectra can be fitted using a simple blackbody and an absorption 
component which has a significant enhancement of neon compared to the solar value. 
We propose that the material in the inter-binary system is significantly enhanced with 
neon. This makes its intrinsic X-ray spectrum virtually identical to RX J0806. 3+1527. 
We re-access the X-ray luminosity of RX J1914. 4+2456 and the implications of these 
results. 

Key words: Stars: binary - general; abundances; individual: - RX J1914. 4+2456, 
RX J0806. 3+1527; X-rays: binaries 



1 INTRODUCTION 

The X-ray source RX J1914.4+2456 (hereafter RX 
J1914+24) has been the subject of much debate since its 
discovery during the ROSAT all-sky survey. A number of 
competing models have been put forward to account for its 
observed properties, but all of them involve a stellar binary 
system. The models can be split into accretion and non- 
accretion models. The non-accretion model is the unipolar 
inductor (UI) model where dissipation of large electrical cur- 
rents heat the magnetic white dwarf (Wu et al 2002). It 
shares many similar observational characteristics to the X- 
ray source RX J0806.3+1527 (hereafter RX J0806+15, see 
Cropper et al 2004a for a review). 

One of the main uncertainties to understanding the na- 
ture of RX J1914+24 is accurately determining its X-ray lu- 
minosity, Lx- In the UI model, Lx is proportional to the rate 
of change of the orbital period and the degree of asynchro- 
nisation between the binary orbital period and the primary 
star (Wu et al 2002, Dall'Osso et al 2007). 

In practise it has been difficult to get an accurate value 
for L x . This is partly due to the fact that RX J1914+24 
is highly reddened. The other factor is that its X-ray spec- 
trum is rather unusual and therefore difficult to determine 
the underlying emission model. X-ray spectra obtained us- 
ing XMM-Newton are not well fitted using a simple absorbed 
blackbody, showing large residuals near 0.7 keV (Ramsay et 
al 2005). Ramsay et al (2005) found that an absorbed black- 
body with a broad emission line centered at 0.59keV gave 



a much improved goodness of fit and for a distance of lkpc 
implied L x ~ 10 35 ergs s _1 . 

Using two further longer series of observations of RX 
J1914+24 also taken using XMM-Newton, Ramsay et al 
(2006) found that an absorbed low-temperature thermal 
plasma model with an edge at 0.83keV gave a significantly 
improved fit compared to the previous best model fit. For 
a distance of lkpc this optically thin emission model gave 
Lx ~ 10 33 ergs s _1 . On the other hand if the distance was 
much lower than 1 kpc (Steeghs et al 2006, Barros et al 2007) 
then Lx could be as low as ~ 3 x 10 31 ergs s _1 - giving a 
range in Lx of 4 orders of magnitude! 

RX J1914+24 has been observed using XMM-Newton 
at 4 separate epochs (Table[l}. An analysis of the data taken 
using the EPIC detectors have been presented in Ramsay et 
al (2005) (from the first two epochs) and in Ramsay et al 
(2006) (from the last two epochs). In this paper we examine 
the data obtained using the RGS detectors; re-examine the 
data obtained using the EPIC detectors using more recent 
calibration data, and also re-examine the models used to fit 
the data. 



2 OBSERVATIONS 

The data were processed using XMM-Newton SAS v7.0 (the 
data presented previousiy were processed using v6.0 and 
v6.5 in Ramsay et al 2005 and Ramsay et al 2006 respec- 
tively). In our analysis we excluded time intervals of high 
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0718 
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6641 


timing 


8933 


sw 


8993 


9880 
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0721 


2003-11-15 


sw 


2888 


timing 


3865 


sw 


3865 


4000 


1.20 x 10" 12 


0880 


2004-09-28 


ff 


11869 


ff 


14954 


ff 


14954 


16954 


1.27 x 10~ 12 


0882 


2004-10-02 


ff 


15274 


ft 


18419 


ff 


18419 


18880 


1.40 x 10~ 12 



Table 1. The log for the observations of RX J1914+24 made using XMM-Newton. We show the mode the detector was in where 'sw' 
refers to 'small window', 'ff' to 'full frame' and the duration of 'good' time — ie excluding time intervals of enhanced background. The 
RGS detectors were configured in 'Spectroscopy' mode. In the last column we show the observed integrated flux in the 0.2-10keV energy 
band using the EPIC pn detector data and fitting an absorbed blackbody plus broad emission line spectral model. 



particle/solar background (a significant issue in the second 
epoch observation). 

For those observations when the EPIC data were in full 
frame mode, we excluded events from the central core of 
the psf (using an aperture of 10 in radius) in order that 
pile-up was not significant. We did not extract spectra from 
the timing mode data since the spectral calibration is not 
as well defined as for the other modes. For the RGS data, 
we extracted spectra which included the source and back- 
ground, and a background spectrum separately. We grouped 
the EPIC spectra so that each bin had a minimum of 40 
counts. Since the RGS spectra from the individual epochs 
were relatively low, we co-added the spectra from the RGS1 
and RGS2 detectors obtained using the 3rd and 4th epoch 
observations (the first two epochs had much shorter expo- 
sures). We used the SAS task rgscombine and then grouped 
each spectrum so that each bin of each spectrum had a min- 
imum of 20 counts per bin. 

To determine the observed X-ray flux at each epoch, 
we used data taken using the EPIC pn detector. We fitted 
the integrated X-ray spectra using an absorbed blackbody 
plus broad emission line. We show the integrated observed 
flux in the 0.2-10keV energy band using in Table Q] (the 
observed flux is only weakly sensitive to the model used). 
This shows that the observed flux varied by 17% between 
the 4 observation epochs. 

The X-ray data folded on the 569 sec period shows a 
distinctive 'on-off ' behaviour, with the X-ray flux being off 
for approximately half the 569 sec period (Cropper et al 
1998). There is a sharp rise in flux which is followed by 
a slower decline from maximum brightness. Ramsay et al 
(2005) showed evidence using the two shorter duration ob- 
servations that the spectrum gets softer during this decline 
phase. Using the 3rd and 4th longer series of observations 
we confirm this finding. Therefore we have obtained a spec- 
trum which covers the 'bright phase' which we define to be 
0=0.08-0.38 where 0=0.0 is defined as the start of the sharp 
rise in X-ray flux. 



3 THE RGS SPECTRA 

We extracted RGS spectra from the bright phase using the 
3rd and 4th epoch observations. In fitting the spectra, we 
used a blackbody, a blackbody plus a Gaussian component 
both in absorption and emission and a multi-temperature 
thermal plasma plus edge model. In the work of Ramsay et 
al (2005, 2006), the absorption model which was used was 
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tbabsvmekaledge 


2.18 


(93) 


tbabsbb 


2.12 


(95) 


tbvarabsbb 


1.28 


(93) 


tbabsbbgau (abs) 


1.20 


(92) 


tbabsbbgau (emi) 


1.14 


(91) 



Table 2. The fits to the bright phase RGS spectrum taken from 
data in XMM-Newton orbits 0880 and 0882. The models noted in 
the first column refer to the models in XSPEC: tbabs - Tubingen 
Boulder absorption ISM model (Wilms et al 2000); tbvarabs - 
Tubingen Boulder absorption ISM model with variable abun- 
dances; bb - blackbody; gau — a Gaussian component in emission 
(emi) and absorption (abs); vmekal - a thermal plasma model 
with non-solar abundances; edge - an absorption edge; In the sec- 
ond column we show the x 2 and (degrees of freedom). 



the 'wa' neutral absorption model found in the XSPEC fitting 
package (Arnaud 1996). Here, we use the Tubingen-Boulder 
absorption ISM model and abundances (Wilms, Allen & Mc- 
Cray 2000) which incorporates advances in atomic cross- 
sections and other physical parameters compared to the 
wa model (Morrison & McCammon 1983). We used this 
model implemented into XSPEC as the tbabs model (which 
assumes an interstellar medium of solar abundance) and the 
tbvarabs model (which allows the abundance of each ele- 
ment to vary). 

We show the goodness of fit to the RGS spectrum using 
the different models in Table [2] As was found by Ramsay et 
al (2005) a simple absorbed blackbody model gives a very 
poor fit. Ramsay et al (2006) found that a low temperature 
thermal plasma plus edge model gave a good fit to the spec- 
trum obtained using the EPIC pn detector. Using the RGS 
data we can rule this model out. The temperature of the 
plasma determined using the EPIC pn detector is very low, 
(< IkeV), which would result in strong X-ray emission lines 
- these lines are not detected in the RGS data. 

This leaves three models - a blackbody with either an 
absorption or emission component, or a blackbody with an 
absorption component which has abundance different to so- 
lar composition. We show the RGS spectrum together with 
the best fit using an absorbed blackbody, where the absorp- 
tion component has variable abundances, in Figure [1] 
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Figure 1. The fits to the RGS spectra bright phase spectra. 
The spectra from each detector (RGS1 - datapoints plotted as 
black circles, RGS2 red points) taken in the 3rd and 4th epoch 
observations. The solid lines show the best fit using an absorbed 
blackbody model (the darker line through the RGS1 data and 
lighter, or red, line for the RGS2 deta), where the absorption 
component has variable abundances. 



4 THE EPIC PN BRIGHT PHASE SPECTRA 

We extracted bright phase spectra from each epoch. We fit- 
ted models consisting of a blackbody and a Gaussian line in 
both emission and absorption, and also a blackbody with an 
absorption component with variable abundances. We show 
the goodness of fits to these spectra using these three mod- 
els in the top panel of Table [3] (For the second epoch ob- 
servation, which had a short good exposure we fixed the 
model parameters at the best fit parameters determined in 
fourth observation apart from the absorption column density 
and the normalisation parameters). We find that the models 
which include a Gaussian in absorption or emission give for- 
mally good fits (95 % confidence) while the model with an 
interstellar component with variable interstellar abundance 
gives a 1-0. 

We show the observed flux in the 0.2-10keV energy band 
and also the unabsorbed flux in the 0.01-10keV energy band 
in the bottom panel of Table [3] The observed flux at each 
epoch as derived from the various models is consistent to 
within 7%. In contrast, the derived unabsorbed flux varies 
by more than an order of magnitude. 

We show the spectral parameters derived for the black- 
body plus a Gaussian component in Table 4. The spectral 
parameters derived using a blackbody plus Gaussian line in 
emission gives consistent results, with a blackbody tempera- 
ture kT ~62 eV and iV H ~ 5 x 10 21 cm -2 . The center of the 
Gaussian component is ~0.65 keV and has a width ~80 eV 
and an equivalent width ~160 eV. In contrast, the spectral 
parameters as derived using the model with the Gaussian 
line in absorption, show a greater spread. However, the tem- 
perature of the blackbody component is hotter than in the 
previous model. The center of the absorption line is ^0.86 
keV and has a width ~0.11 keV. 

Since the variable abundance absorption model was not 
used by Ramsay et al (2005, 2006) we now go onto inves- 
tigate this model in more detail. We fit the spectra using 



an absorbed blackbody model where we allowed the abun- 
dance of the absorbing material to vary from solar. We did 
this by allowing each element to vary from solar one at a 
time. If there was no significant change in the goodness of 
fit, we re-set the abundance for that element back to solar 
and fixed this parameter. We did this for the EPIC pn and 
RGS spectra from XMM-Newton orbits 0880 and 0882 indi- 
vidually. We show the spectral parameters for these fits in 
Table [S] We find that their X-ray spectra can be well fit- 
ted using blackbody with an interstellar absorption model 
which has significantly increased amounts of Neon. There is 
some evidence that Iron could have non-solar abundances. 

We also fitted the EPIC pn spectra taken from XMM- 
Newton orbits 0718 and 0721 using this same model. As 
mentioned before, for the data taken in the orbit 0721 ob- 
servations, we fixed the spectral parameters at their best 
fit parameters determined in orbit 0882 (apart from the ab- 
sorption column density and the blackbody normalisation). 
We show the spectral parameters derived from these spec- 
tra in Table [5] We find these spectra are also consistent with 
the abundance of neon being significantly enhanced in these 
epochs as well. 



5 AN EVALUATION OF THE SPECTRAL 
MODELS 

We know go onto to discuss the physical plausibility of the 
three spectral models which we have used. Namely, the ab- 
sorbed blackbody with an absorption component and an 
emission component, and the blackbody with absorption 
component with non-solar abundances. 

5.1 A blackbody with Gaussian absorption 
component 

Isolated neutron stars (INS) with high magnetic fields have 
been found to show absorption lines which have been at- 
tributed to either a proton cyclotron line or an electron 
cyclotron line, see Zane et al (2005), Haberl (2007) and 
Schwope et al (2007). The similarity between the spectral 
parameters of RX J1914+24 and some INS is quite striking. 
For instance, the width of the line and its equivalent width 
as measured in RX J1914+24 is very similar to that of RBS 
1223 (Schwope et al 2007). Also the variation in the observed 
flux in the energy range 0.35-1. 5keV of the INS RX J0720.4- 
3125 is 20% (cf 17 % between the XMM-Newton observa- 
tions of RX J1914+24) showing that non-accreting sources 
can show a significant variation in their X-ray flux (the data 
on RX J0720.4-3125 were extracted from the XMM-Newton 
archive) . 

There are, however, some very significant differences be- 
tween the observed properties of RX J1914+24 and that of 
INS. The first is that their spin-periods are in the range 
~ 3 — 12 sec - these are very much shorter than the period 
seen in RX J1914+24 (569 sec). The second is the luminos- 
ity difference - INS show X-ray luminosities J;10 31 ergs s . 
Using the inferred unabsorbed fluxes as derived using the 
blackbody plus Gaussian absorption line, we find that RX 
J1914+24 would have to be at a distance of 20 pc to give 
a comparable luminosity. If RX J1914+24 was so close, we 
would expect to detect a significant proper motion which 
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Table 3. Top panel: The fits to the bright phase EPIC pn spectrum obtained using XMM-Newton at four epochs. The models noted 
in the first column are defined in the caption for Tabic 2. In the following columns we show the goodness of fit (x 2 and degrees of 
freedom). Bottom panel: we show the observed flux, Flux°, in the 0.2-10keV energy band, and the implied unabsorbed flux, Flux 11 , in 
the 0.01-10keV energy band for all four epochs and the three spectral models. 



XMM N-n kT bb E a EW 

Orbit (10 21 cm- 2 ) (eV) (keV) (kcV) (eV) 



Gaussian component in emission 

0718 4.2+0; 5 , 64.9+1; 9 0.64+°;° 2 0.0751°;°}* 182+^' 

0880 5 - 7 -o.2 57 - 5 l 2 '.o 0.63±0.01 0.093tomt 152 -a 

0882 S-ito.T- 63.9^47 0.68±0.01 0m±g°l 147l| 

Gaussian component in absorption 

0718 2.9+°"3 110+18 086+0^ o.lll±°-gg| -216lg, 

0880 4.9l^ 80.7115-5 0.87l» : ^ 0.12±°-°J -148lg* 

0882 6.2l°;| 70.6l|J 0.84l°;° 2 0.1051°;°^ -W4±H 



Table 4. The spectral parameters for fits to the bright phase spectra using an absorbed blackbody model where we have allowed the 
abundance of the absorption component to vary from solar. The duration of the EPIC pn spectrum from XMM-Newton orbit 0721 was 
relatively short and hence the spectral parameters are not strongly constrained. 



has not been observed (Israel et al 2002). The third differ- 
ence is the change in the period. The 569 sec period in RX 
J1914+24 has been found to be decreasing, while the period 
of the two INS which have been found to show a change 
in their period are increasing (Cropper et al 2004b, Kaplan 
& van Kerkwijk 2005a, Kaplan & van Kerkwijk 2005b). A 
fourth difference is the optical brightness, with INS being 
typically B ~26 (see the compilation in Haberl 2007), while 
RX J1914+24 is B ~21. We conclude RX J1914+24 is not 
an isolated neutron star. At this point, whilst we cannot rule 
out the presence of absorption features in the X-ray spec- 
trum of RX J1914+24, we do not consider it the most likely 
model to explain its X-ray spectrum. 

5.2 A blackbody with Gaussian emission 
component 

A blackbody with an additional broad emission line does, on 
first sight, seem rather contrived. However, such a feature 
has been claimed to be present in the relatively low reso- 
lution ASCA spectra of a number of X-ray ultra-compact 
binaries with neutron star primaries (eg Juett, Psaltis & 
Chakrabarty 2001). The line center of the emission line is 
remarkably similar in RX J1914+24 and the four sources 



shown described in Juett et al (2001). One notable differ- 
ence is that the temperature of the blackbody component 
in the X-ray UCBs is much hotter - several hundred eV as 
opposed to ~60 eV - and the fact that they are hard X-ray 
sources being detected at energies up to many 10's of keV. 
This is the result of the primary being a neutron star as 
opposed to a white dwarf. 

It was claimed that this broad emission feature was 
due to the superposition of a number of unresolved emis- 
sion lines. However, when one of the sources observed using 
ASCA was observed using the Chandra Low Energy Grating 
Spectrometer, it failed to detect any emission features which 
could give rise to a broad emission feature in low resolution 
spectra. We now address one possible reason for this. 

5.3 A blackbody with variable abundances in the 
absorption model 

Juett et al (2001) found that good fits to Chandra spectra 
of neutron star UCBs were obtained if the absorption model 
had non-solar abundances. In particular, they found a high 
relative abundance of neon and suggested that this over- 
abundance was located in the inter-binary system. However, 
further work (eg Juett & Chakrabarty 2005) showed that for 
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Table 5. The spectral parameters for fits to the bright phase spectra using an absorbed blackbody model where we have allowed the 
abundance of the absorption component to vary from solar. Since the duration of the spectrum taken in orbit 0721 was short, the spectral 
parameters were fixed at their values determined in orbit 0882 apart from the total absorption column density and normalisation of the 
blackbody component. 



individual sources, the Ne/O ratio showed evidence for vari- 
ability from epoch to epoch, which they attributed to source 
variability. This implied that the abundances could not be 
used to determine the composition of the mass donating star. 

There is a clear similarity between the neutron star 
UCBs described by Juett et al and RX J1914+24. In each 
observation of RX J1914+24, there is clear evidence that 
the absorption component has an overabundance of neon. 



6 DISCUSSION 

For the reasons outlined in £15.11 we rule out RX J1914+24 
being an isolated neutron star. Since all the known neutron 
star UCBs have X-ray emission extending up to many 10's of 
keV we also rule out an accreting neutron star UCB model. 
We cannot rule out that a neutron star is in a binary system 
where a secondary star was not filling its Roche Lobe. In this 
scenario, an X-ray bright system would have to be powered 
by UI. 

It is highly unlikely that the line of sight absorption to 
RX J1914+24 has a chance enhancement of neon. It is much 
more likely that this overabundance is concentrated in the 
binary system. Juett & Chakrabarty (2005) noted that for 
some neutron star UCBs the Ne/O abundance varied from 
epoch to epoch and hence the observations could not be used 
to determine the abundance of the secondary, mass-donating 
star, in the binary system. In the case of RX J1914+24 there 
is clear evidence for a significant over-abundance of neon in 
the absorption component at each epoch. At this stage it 
is not clear if this over-abundance is due to circumbinary 
material left over from a previous stage in the binary forma- 
tion process or can give us a direct insight into the chemical 
composition of the secondary star (if accretion is occurring). 

What are the implications of our findings regarding the 
X-ray luminosity of RX J1914+24? Steeghs et al (2006) dis- 
cuss the extinction and distance estimates to RX J1914+24. 
While the distance is rather uncertain, it is likely that it 
is greater than ~1 kpc. We can rule out the lower esti- 
mates (Lx ~ 10 33 ergs s~ 1 for a distance of 1 kpc) which 
were derived using a low temperature thermal plasma model. 
Taking the unabsorbed bolometric fluxes derived using the 
blackbody with absorption component with variable abun- 



dances and assuming a distance of 1 kpc we find Lx = 
2 x 10 34 - 1.6 x 10 35 ergs s _1 . 

Dall'Osso et al (2007) made a detailed investigation 
of the UI model in the context of RX J1914+24 and 
RX J0806+15. They predicted that for low luminosities, 
Lx ~ 10 33 ergs s" 1 , the asynchronism between the orbit 
and the magnetic star in RX J1914+24 would have to be 
q ~ 0.9—0.98, where a — uti/(j , and u>i is the rotation fre- 
quency of the primary star and uj is the orbital frequency. 
Unless RX J1914+24 was located at a distance significantly 
less than lkpc, we can rule these low values of asynchronism. 
For high luminosities (L x = 10 34-35 ergs s _1 ), Dall'Osso et 
al (2007) predicted that an asynchronism of a few was re- 
quired (a ~ 4). For an observed period of 569 sec, a = 2 — 10 
gives a predicted period of ~60-300 sec. There is no evidence 
for power at these periods in the power spectra of the X-ray 
light curves (Ramsay et al 2006). 



7 SUMMARY 

Until now the nature of the emission source that powers the 
X-ray spectrum of RX J1914+24 has been far from clear. 
In this paper we have shown that it can be well modelled 
using a simple blackbody model with an absorption com- 
ponent which has non-solar abundances, in particular, an 
enhancement of neon. 

Since the X-ray light curves of RX J1914+24 and RX 
J0806+15 are practically identical, it suggests that their X- 
ray emission source is the same. The fact that their X-ray 
spectra were apparently different (albeit both being soft) 
was therefore perplexing. Our result showing that the emis- 
sion source is the same for both RX J1914+24 and RX 
J0806+15 is therefore very attractive. Indeed their temper- 
atures are virtually identical - we obtain a mean value of 
kT ~ 67 eV for RX J1914+24 compared to kT ~ 65 eV for 
RX J0806+15 (Israel et al 2003). 

The difference between the X-ray spectrum of RX 
J1914+24 and RX J0806+15 is that the absorption compo- 
nent of RX J1914+24 has enhanced neon abundance. A fur- 
ther investigation of the optical spectrum of RX J1914+24 
to search for neon features is strongly encouraged. 
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